Astronomy & Astrophysics manuscript no. ngc3516-rgs 


©ESO 2010 


February 18, 2010 





The warm absorber and X-ray variability of the Seyfert 1 galaxy 
NGC 3516 as seen by the XMM-Newton RGS 

M. Mehdipour* , G. Branduardi-Raymont, and M. J. Page 



o 

00 



o 

U 

6 



> 

O 

(N 

m 

(N 
O 
O 



X 



Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, United Kingdom 
Received August 2009 / Accepted February 2010 

ABSTRACT 

Aims. We present a new analysis of the soft and medium energy X-ray spectrum of the Seyfert 1 galaxy NGC 3516 taken with the 
Reflection Grating Spectrometer (RGS) and European Photon Imaging Camera (EPIC) on board the XMM-Newton observatory. We 
examine four observations made in October 2006. We investigate whether the observed variability is due to absorption by the warm 
absorber and/or is intrinsic to the source emission. 

Methods. We analyse in detail the EPIC-pn and RGS spectra of each observation separately. 

Results. The warm absorber in NGC 3516 is found to consist of three phases of ionisation, two of which have outflow velocities of 
more than 1000 km s . The third phase (the least ionised one) is much slower at 100 km s . One of the high ionisation phases, with 
log £ of 2.4, is found to have a partial covering fraction of about 60%. It has previously been suggested that the passage of a cloud, part 
of a disk wind, in front of the source (producing a change in the covering fraction) was the cause of a significant dip in the lightcurve 
during one of the observations. From our modelling of the EPIC-pn and RGS spectra, we find that variation in the covering fraction 
cannot be solely responsible for this. We show that intrinsic change in the source continuum plays a much more significant role in 
explaining the observed flux and spectral variability than originally thought. 

Key words, galaxies: active - galaxies: Seyfert - galaxies: individual: NGC 3516 - X-rays: galaxies - techniques: spectroscopic 



1. Introduction 

Recent X-ray observations of many Seyfert 1 galaxies and some 
quasars show evidence of absorbing phot oionised gas in our 



line of sight towards the active nucleus (e.g. Kaastra et alJ 



Reeves et af1| 2003b iBlustin et all 120051: ICostantini etaTI 



2000; 



2007 



Kaastra 2008). This gas is normally found to be outflowing 



with often several distinct components, and is referred to as a 
"warm absorber"; its absorption signatures are clearly seen in 
the high resolution X-ray spectra. AGN outflows, in general, are 
believed to have important astrophysical implications. For ex- 
ample, strong outflows from the AGN can affect the growth of 
the supermassive black hole at its core and enrichment of the 
ISM/IGM could alter the course of the host galaxy evolution 
through feedback processes. The extent of the outflow contri- 
bution depends on its kinetic luminosity and the mass outflow 
rate. To reconstruct the kinetic luminosity of the outflow, the 
ionisation state and structure of the warm absorber need be to 
determined from photoionisation modelling of the high resolu- 
tion X-ray spectra. 

The first w arm absorber identification was reported by 
iHalpernl dl984 from an observation of the QSO MR2251-178 
with the Einstein Observatory. Then in the days of ROSAT and 
ASCA, warm absorbers were found through the detection of 
what was interpreted as broad and blended continuum absorption 
edges of ions such as O vn and O vm. For exa mples of early stud- 
ies of AGN warm absorbers using ROS AT seelNandra & Pounds! 
(1 19921). iNandra et all (fl99l and jTurner et al.1 dl993l): f o r thos e 
using ASCA, see Reynolds] d 19971) and iGeorge et al.1 d!998l) . 
However, the instruments used at that time had insufficient spec- 



Send offprint requests to: M. Mehdipour 
* e-mail: missagh . mehdipour@ucl .ac.uk 



tral resolution to measure any details, such as the presence 
and parameters of narrow absorption lines. With the advent of 
high resolution X-ray spectrometers onboard XMM-Newton and 
Chandra, the availability of high-resolution spectra has allowed 
a major leap forward in the study of AGN outflows in the last 
decade, starting with the original discovery of numerous blue- 
shifted absorp tion lines from pho toionised, outflowing gas in 
NGC 5548 by iKaastra et all dioOOl) . the subsequent identifica- 
tion of similar features in a number of AGN and the discovery 
of an unresolved transiti on array (UTA) of l ow ionised iron ions 
in IRAS 13349+2438 bv lBehar et al1d200lh . 

NGC 3516 is a Seyfert 1.5 SB0 galaxy at a redshift of 
0.008836 dKeelll 19961) in the constellation of Ursa Major. Even 
from early X-ray observations of this object, signatures of a 
multi-phase warm absorber have been evident. From simultane- 
ous far-UV and ASCA X-ray observations in 1995, a warm ab- 
sorber with at lea st two absorption components was reported by 
Kriss et al.l d!996l) . The two components differ by a factor of 8 in 
ionisation parameter; the more highly ionised has a column den- 
sity tw ice as large as the less ionised component. iNetzer et al.l 
d2002l) investigated spectral variations of NGC 3516 over a pe- 
riod of 7 years, by using archival ASCA and early Chandra ob- 
servations. They reported a large drop in flux (factor of ~ 50 at 
1 keV) between an ASCA observation i n 1994 and the Chandra 
observation in 2000. lNetzer et al.l(l2002l) concluded that the vari- 
ations in the observed flux and spectra at these epochs were con- 
sistent with a constant column density of line-of-sight material 
reacting to changes in the ionising continuum. 

NGC 3516 was observed twice by XMM-Newton in April 
and November 2001; both observati ons were pa r tially over- 
lapping with Chandra observations. iTurner et al.l ((2002) pre- 
sented results from the simultaneous Chandra High Energy 
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Transmission Grating (HETG) and XMM-Newton observations 
made in November 2001: analysis of the Fe Ka regime showed 
evidence of several narrow emission features and of rapid evo- 
lution of the F e Ka line during t he observation. From the 2001 
observations. iTurner et al.l (12005) reported the presence of three 
distinct zones (phases) of gas with different ionisation parame- 
ters covering the active nucleus. The ionisation parameter (£) is 
defined as 

I*» (1) 
nr L 

where Li on is the luminosity of the ionising source over the 1— 
1000 Ryd band (in erg s" 1 ), n the hydrogen number density (in 
cm -3 ) and r the distance between the ionised gas and the ionising 
source (in cm). Therefore, £ is in units of erg cm s" 1 . 

The warm absorber phases found in lTurner et alj d2005l) are: 
a low-ionisation UV/X-ray absorber with log £ 0.5 and hy- 
drogen column density of A^h ~ 0.5 x 10 22 cirT 2 ; a more highly 
ionised gas with log ^ ~ 3.0 and Nh ~ 2 x 10 22 crrT 2 outflowing 
at a velocity of ~ 1 100 km s" 1 ; and a phase with log £ ~ 2.5 and 
a very large hydrogen column de nsity of Nn ~ 25 X 10 22 cm" 2 
covering ~ 50% of the continuum. Turner et al.l (|2005 ) found the 
spectral variability in the 2001 observations to be consistent with 
the ionisation-state of the absorbing gas layers responding to the 
continuum flux variation. 

More recently, four observations of NGC 3516 were per- 
formed by XMM-Newton in October 2006, interwoven with 
Chandra observations. The analysi s of the EPIC-p n and 
Chandra HETG data is published in ITurner et al.l 12008). They 
discovered a previously unknown ionisation phase in the Fe Ka 
regime, with log £ ~ 4.3 and A^h ~ 26 x 10 22 cm" 2 . In line 
with the results from the 2001 observations, the presence of a 
phase with a covering fraction of ~ 50% was confirmed. This 
phase was found to have log £ ~ 2.2 and A^h ~ 20 x 10 22 cm" 2 . 
The source showed significant flux variability between observa- 
tions, especially a d ip in the lightcurve of the third observation. 
ITurner et al.l ([2008) concluded that changes in the covering frac- 
tion of the phase partially covering the continuum provide a sim- 
ple explanation of the dip in the lightcurve. They interpreted this 
as an eclipse of the continuum due to passage of a cloud across 
the line of sight over half a day. 

In this paper, we present a new analysis of the EPIC-pn, 
RGS and OM data from the four XMM-Newton observations 
of October 2006, with the emphasis put on the RGS spectra. 
Section 2 describes the observations and data analysis and Sect. 
3 focuses on the lightcurves; the spectral modelling is described 
in detail in Sects. 4 and 5; we discuss our findings in Sect. 6 and 
give concluding remarks in Sect. 7. 

2. Observations and data analysis 

We extracted the EPIC-pn, RGS and OM data for the four 
2006 observations from the XMM-Newton Science Archive 
(XSA). The XMM-Newton observation IDs are 04012010401 
(51.6 ks), 04012010501 (68.8 ks), 04012010601 (68.2 ks) and 
04012011001 (67.8 ks). Hereafter, we refer to them as Obs. 1, 
Obs. 2, Obs. 3 and Obs. 4, respectively. 

All the Observation Data File (ODF) data were processed us- 
ing SAS v7. 12, unless specified otherwise in the text. For EPIC- 
pn, RGS 1 and RGS2, periods of high flaring background were 
filtered out before extracting scientific products. For the EPIC- 
pn, this was done by creating Good Time Intervals (GTIs), and 
using them in conjunction with the internal GTI tables to re- 
move periods with single event count rates exceeding 0.4 counts 



per second in the energy range E > 10 keV. For the RGS instru- 
ment, in addition to the internal GTIs, we constructed an extra 
GTI table to filter out periods with count rates above 0.4 counts 
per second in CCD number 9, which is the closest to the optical 
axis of the telescope and therefore the most affected by back- 
ground flares. 

The spectra and lightcurves of EPIC-pn, operating in the 
Small Window mode, were extracted from a circular region of 
40" radius centred on the source. The background was extracted 
from a source-free region of the same radius, located at the same 
RAWY as the source on the same chip. The EPIC-pn showed no 
evidence of pile-up. The EPIC-pn lightcurves were background- 
subtracted and corrected (using the epiclccorr SAS v8 task) 
for various effects on the detection efficiency such as vignetting, 
bad pixels, chip gaps, Point Spread Function (PSF) variation 
and quantum efficiency; the task also makes corrections which 
vary with time affecting the stability of the detection, such as 
dead time and GTIs. The RGS instruments were operated in the 
standard Spectro+Q mode. The first and second order spectra of 
RGS 1 and RGS2 were extracted, and the response matrices were 
generated using the rgsproc processing task. 

The OM was operated in Image+Fast mode in all observa- 
tions, as well as in Science User Defined mode in the last two ex- 
posures of Obs. 1, Obs. 3 and Obs. 4. The images in all four ob- 
servations were taken using the U filter (3000-3900 A). The OM 
Image mode data were processed using the omi chain pipeline. 
The resulting images are each of 1400 s exposure. The FWHM 
of the source is about 1.6" which is consistent with the OM 
on-board PSF FWHM of 1.55" in the U filter. We performed 
photometry on each image in a fully interactive way using the 
omsource program: we carefully selected the source and back- 
ground regions to extract the count rate, and applied all the nec- 
essary corrections such as for the PSF and coincidence losses, in- 
cluding time-dependent sensitivity (TDS) corrections. The OM 
lightcurves were extracted from a circle of 8 pixels radius (1 
pixel = 0.48") centred on the source nucleus. The background 
was extracted from a source-free region of the same radius. 

All the results presented in this work are from simultaneous 
fits of the EPIC-pn and the first-order RGS1 and RGS2 spec- 
tra, unless specified otherwise in the text. The EPIC-pn spectra 
are in the range 0.2-10 keV and are binned to a minimum of 200 
counts per bin. The RGS data are in the range 6-38 A (0.33-2.07 
keV), binned by a factor of 3 (this corresponds to an average bin 
size of 0.04 A, which still over-samples the RGS resolution ele- 
ment of ~ 0.07 A FWHM). All the spectra shown in this paper 
are background-subtracted and are in the observed frame. The 
spectral analysis and mo delling was done us ing the SPEX (ver- 
sion 2.01.02) fitting code dKaastra et all 1996). All the parameter 
errors quoted in this paper correspond to a Ax 2 of 2. 



3. X-ray and UV lightcurves 

Figure Q] shows the OM (UV, top) and EPIC-pn (X-ray) 
lightcurves of the four observations. The source shows signif- 
icant X-ray flux variability between observations, notably the 
flare in Obs. 1 and the dip in Obs. 3. From visual inspection 
of the X-ray lightcurves, the variations in the 0.2-2.0 keV (soft 
X-ray) and 2.0-10.0 keV (hard X-ray) bands follow each other 
closely; however in Obs. 3 the ratio of soft to hard X-ray count 
rates clearly increases as the observation progresses indicating 
spectral variability; some evidence for a softer spectrum when 
brighter is also present in Obs. 1 and Obs. 4. 
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Fig. 1. Upper panel: Top lightcurves: the OM background- 
subtracted U band (3000-3900 A) lightcurves of Obs. 1 (red), 
Obs. 2 (purple), Obs. 3 (black) and Obs. 4 (blue); each count 
rate is scaled down by a factor of 4 for clarity of display and cor- 
responds to an OM exposure of 1400 s. Middle lightcurves: the 
EPIC-pn background-subtracted lightcurves of the four observa- 
tions in 200 s time bins over 0.2-2.0 keV. Bottom lightcurves: 
the EPIC-pn background-subtracted lightcurves of the four ob- 
servations in 200 s time bins over 2.0-10.0 keV. Lower panel 
(green): the ratio of 0.2-2.0 keV (soft X-ray) to 2.0-10.0 keV 
(hard X-ray) lightcurves. 

In contrast to the X-ray, the UV lightcurve shows no signifi- 
cant flaring in Obs. 1 and there is not a significant dip in the UV 
flux in Obs. 3; it is not unus ual to see variability in the X-rays, 
and none in the UV (e.g. see iBlustin et al.l l2~002 for analysis of 
NGC 3783). 

4. Preliminary modelling of the X-ray continuum 
and Fe Ka line 

Figure |2]depicts an overview of the EPIC-pn and RGS spectra of 
the 4 observations: while Obs. 2 and 4 display very similar spec- 
tra (so much so that they are practically indistinguishable), that 
of Obs. 1 lies slightly above them, and that of Obs. 3 is clearly 
much fainter than all the others. The features of the spectra, in 
particular during the dip in Obs. 3, are discussed in the following 
sections. 

We started modelling the 4.0-10.0 keV EPIC-pn spectrum 
of Obs. 2 (which is the longest of the 4 observations) by a sim- 
ple power-law (pow) redshifted to z — 0.008836. The transmis- 
sion of the Galactic neutral absorption was included by apply- 



[__, — , — , — , — — i — , — i ■ ■ i — , — , — , — — , — , — , — , — i—i — , — , — , — i — , — i — , — , — i — , — , — I 

5 10 15 20 25 30 35 

Observed Wavelength (A) 

Fig. 2. The EPIC-pn spectra of Obs. 1 (red), Obs. 2 (purple), 
Obs. 3 (black) and Obs. 4 (blue) binned to a minimum of 200 
counts per bin. In the inset are the fluxed, combined RGS 1 and 
RGS2 spectra binned by a factor of 3. Note that Obs. 2 and Obs. 
4 are practically indistinguishable. 
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Observed Energy (keV) 

Fig. 3. The Obs. 2 EPIC-pn power-law fit (including the 
Gaussian Fe Ka line and Galactic absorption) over the 4.0-10.0 
keV range. The fit is extrapolated to lower energies, displaying 
the presence of a soft excess and additional absorption in the 
0.6-2.5 keV range. 



ing the hot model (collisional ionisation equilibrium) in SPEX. 
The Galactic H i column density in our line o f sight was fixed to 
N H = 3.45 x 10 20 cnr 2 (iKalberla et al.ll2005l) and the gas tem - 
perature to 0.5 eV to mimic a neutral gas dKaastra et al.l l2009). 
The prominent Fe Ka emission line at ~ 6.4 keV was mod- 
elled by a simple Gaussian profile (gaus) with its width left 
free to vary. The parameters obtained from the preliminary mod- 
elling of the continuum and Fe Ka line are shown in Table Q] 
Interestingly, we find that whereas the power-law flux falls sig- 
nificantly in Obs. 3, that of the Fe Ka emission line remains 
unchanged compared to the other observations. 

Figure |3] shows the 4.0-10.0 keV EPIC-pn fit, extrapolated 
to lower energies, displaying the presence of a soft excess and 
additional absorption in the 0.6-2.5 keV band. The latter is most 
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Table 1. Power-law and Fe Ka emission line parameters, obtained from preliminary EPIC-pn fits (including the Galactic absorption) 
of all 4 observations over the 4.0-10 keV energy range. 





Power-law: 










FeKor: 










Obs. 


Photon Index (r) 


Normalisation " 




E 


C 


Flow v d 


FWHM " 


O", 


/ 


Normalisation 8 


^/ d.o.f. 


1 


1.85 ±0.02 


2.9 ±0.1 


6.39 


6.37 ± 


0.02 


+900+™° 


280 ± 70 


6000 ± 


2000 


1.1 ±0.1 


1.54/330 


2 


1.82 ±0.02 


2.5 ±0.1 


6.39 


6.38 ± 


0.01 


+500^° 


290 ± 40 


6000 ± 


1000 


1.3 ±0.1 


1.42/361 


3 


1.75 ±0.02 


1.9 ± 0.1 


6.39 


6.38 ± 


0.01 


+500+ 400 


290 ± 40 


6000 ± 


1000 


1.3 ±0.1 


1.58/324 


4 


1.84 ±0.02 


2.6 ±0.1 


6.39 


6.40 ± 


0.01 


-500+ 500 


260 ± 40 


5000 ± 


1000 


1.2 ±0.1 


1.38/311 



" 10 S1 photons s -1 keV~' at 1 keV. 

* Theoretical rest frame wavelength in keV. 

Wavelength in the rest frame of NGC 35 16 in keV. 

d kms" 1 . 

' eV. 

f RMS velocity in km s" 1 (cr r = FWHM/ Vhi256). 

* 10 49 photons s~'. 



likely due to the warm absorber and will be modelled in the next 
section, by fitting the RGS and EPIC-pn data simultaneously. 

To include modelling of the soft excess, we refitted the EPIC- 
pn spectrum of Obs. 2 over the 0.2-10.0 keV energy range. We 
found that the power-law continuum model cannot fit the spec- 
trum at all; a reduced Chi-squared (xl) value of 143 is obtained 
for 1111 degrees of freedom (d.o.f.). We then added a modified 
black body component (mbb) to the continuum in addition to the 
power-law. The mbb model describes the spectrum of a black 
body modified by coherent Compton scattering, a nd is often used 
to rep resent the accretion disk spectra of AGN dKaastra & Bard 
1989). The addition of the mbb component (kT = 128 eV) 
caused the xl to fall from 143 (1 1 1 1 d.o.f.) to 21.0 (1 109 d.o.f.). 
This is not yet a satisfactory fit, which shows the need for fur- 
ther modelling of the absorption in the 0.6-2.5 keV range (i.e. 
the warm absorber); this is discussed in the following section. 

5. Spectral modelling of the warm absorber 

In the following subsections we describe the process of fitting 
the spectra with a model including warm absorption. We start 
by fitting in detail the spectrum of Obs. 2, which represents an 
"average" state of the source. In Sect. 15.11 we fit only the RGS 
spectrum of Obs. 2 and in Sect. 15. 21 both the RGS and EPIC-pn 
simultaneously; we then proceed to fit the spectra of all the other 
observations. 

5. 1 . Spectral fit using a thin slab absorption model 

We began the analysis of the warm absorber by identifying and 
characterising the main absorption lines in the RGS spectrum 
of Obs. 2, shown in Fig. [4] The slab model in SPEX calculates 
the transmission through a slab of material with adjustable ionic 
column densities. We started by adding a slab component to the 
best fit continuum model obtained from the EPIC-pn. All the 
continuum parameters except the power-law photon index were 
left free. The photon index was let free only after an overall good 
fit was obtained. This helps prevent being caught in false xl min- 
ima and obtaining a photon index incompatible with the EPIC- 
pn during the^J minimisation. This method of fixing the power- 
la w photon index until clos e to the end of the fitting is also used 
in Steenbrugge et al.l (l2003h . 

Initially all the ions in the slab component have negligible 
column densities. We fitted the column density of each ion one 




Observed Wavelength (A) 

Fig. 4. The RGS spectrum of Obs. 2, fitted using the slab model 
described in Sect. I5.ll Some of the strongest absorption features 
are labelled. 

Table 2. Flow and RMS (broadening) velocities of the ions with 
the strongest absorption lines in the Obs. 2 RGS spectrum, de- 
termined from the slab model best fit. 



Ion 


Flow v 


a 


RMS 


v" 


Observed Transitions 


Cvi 


-800d 


: 100 


300 ± 


100 


Ly-a-, Ly-/ 


?, Ly-y 


Nvi 


-700i 


:300 


600 ± 


200 


He-or, He- ; 


8, He-y 


N vn 


-800 ± 


:300 


900 ± 


300 


Ly-a-, Ly-/ 


1, Ly-y 


Ovu 


-900 ± 


100 


300 ± 


100 


He-or, He- ; 


8, He-y, He-c) 


O viii 


-1000 


± 100 


370 ± 


40 


Ly-a-, Ly-/ 


?, Ly-y, Ly-S 


km s . 



by one, starting from the ions responsible for the strongest ab- 
sorption lines in the spectrum. All the ions in a slab component 
have the same flow and RMS (broadening) velocities. So after 
fitting the column densities of all the ions, we decoupled the 
velocities of some of the ions by creating a separate slab com- 
ponent for each ion. Table|2]shows the velocities determined for 
these ions. Only decoupling the velocities of the ions with strong 
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absorption lines improves the fit. Decoupling velocities of the 
ions with weak lines does not improve the fit because large er- 
rors are associated with these ions. The ions with weak lines, that 
do not appear in Table [2] are modelled as one slab component 
with an outflow velocity of 900 ±100 km s _1 and RMS velocity 
of 900 ± 100 km s -1 . 

We only used the slab model to get some understanding of 
the scale of velocity and column density of the ions, especially 
those with the clearest absorption lines in the RGS spectrum (e.g. 
O vm). To analyse the warm absorber in detail, we used a more 
realistic model described in the following subsection. 

5.2. Spectral fit using a photoionised absorption model 

The xabs model in SPEX calculates the transmission through 
a slab of material, where all ionic column densities are linked 
in a physically consistent fashion through a grid of XSTAR 
dBautista & Ka llman 2001 ) photoionisation models. We used the 
default xabs model in which the Spectral Energy Distribution 
(SEP) is based o n the one from NGC 5548 as used in 
Steenbrug ge et al.1 d2005l) . This is a good approximation for 
NGC 3516 since both Seyfert 1.5 AGN have very similar SEDs 
as shown in NASA/IPAC Extragalactic Database (NED). The ad- 
vantage of the xabs model over the slab model is that all relevant 
ions are taken into account, including those, which would be de- 
tected only with marginal significa nce or not detected at all using 
the slab model dKaastra et al.ll2009l) . 

We started with the EPIC-pn best fit continuum and Fe Ka 
model established in Sect. |4] leaving the power-law, mbb and Fe 
Ka parameters free (except the power-law photon index, which 
was kept fixed to the value given in Table Q] until close to the 
end of the fitting procedure). We applied the xabs model to the 
continuum to fit the ionisation parameter (£), the hydrogen col- 
umn density (Nn), the flow and RMS velocities in the EPIC-pn 
and RGS spectra of Obs. 2 simultaneously. The e lemental abun - 
dances were fixed at the proto-solar values of lLodderi ([2003). 
Because we were fitting EPIC-pn, RGS1 and RGS2 spectra 
simultaneously, we allowed the relative normalisations of the 
EPIC-pn and the two RGS instruments to be free parameters. We 
fixed the instrumental normalisations of RGS1 to 1, and freed 
the normalisation of the RGS2 and EPIC-pn instruments. By do- 
ing this, problems with differences in rel ative normalisation s be- 
tween different instruments are avoided (Blu stin et al.f 2003). 

We then included one xabs component in the model. This 
corresponds to phase B in Table [4] where the phases are sorted 
by increasing value of Adding the first xabs component im- 
proved xl from 25.3 (2801 d.o.f.) to 3.91 (2797 d.o.f.), by fit- 
ting Nn, flow and RMS velocities. Our model still required 
more absorption at ~ 11-15 A (expected to be due to hot iron, 
Fe xvn-Fe xxiv) and ~ 16— 17 A (most likely due to cold iron, 
Fei-Fexiv). So to produce a better fit to the data, we intro- 
duced a second xabs component to our model. This corresponds 
to phase C in Table [4] The addition of phase C, which is the 
one with the highest ionisation parameter, further improved xl 
from 3.91 (2797 d.o.f.) to 1.88 (2793 d.o.f.). We found that af- 
ter adding two xabs components, absorption by iron, especially 
the M-shell iron (Fe i-Fe xiv) forming the Unresolved Transition 
Array (UTA) at ~ 16-17 A dBehar et al.l200ll) . was still not prop- 
erly modelled. To model the UTA, a phase with ionisation lower 
than the two previous ones identified was required. Therefore, a 
third xabs component (phase A) was added to our model. The 
addition of phase A improved^ from 1.88 (2793 d.o.f.) to 1.63 
(2789 d.o.f.) as the UTA was fitted well. At this point, the power- 
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Observed Wavelength (A) 

Fig. 5. The Obs. 2 best fit xabs model, fitted simultaneously to 
the RGS and EPIC-pn spectra. The RGS data are shown in blue, 
the EPIC-pn data in red and the models in black. Close-ups of the 
RGS spectrum, showing the spectral features more clearly, are 
presented in Fig. [6] Note that both RGS and EPIC-pn spectra in 
this figure have been divided by the instrument's effective area, 
thus are plotted on the same scale. 

law photon index (F) was left free to vary in the fitting procedure 
for Obs. 2, but neither the value of T (1.82) nor^ changed sig- 
nificantly. 

To attempt to improve further the quality of our fit, we freed 
the abundances of only the elements with the strongest absorp- 
tion features in the RGS spectrum: N, O and Fe. In order not to 
increase the number of free parameters unnecessarily, we cou- 
pled the abundances of all three phases. The precise abundances 
of elements in each phase cannot be accurately determined since 
letting them all free results in some cases in unphysical values. 
Assuming that different phases have the same abundances is al- 
ready a more progressive approach than the alternative of assum- 
ing cosmic abundances for all the phases. The abu ndance cou- 
pling approach is also used by iKaastra et al.l (120031) where they 
fixed the abundances of a cool gas component to those of a hot 
one, in order to reduce the number of free parameters. Freeing 
the N, O and Fe abundances, produced a better fit as xl fell from 
1.63 (2788 d.o.f.) to 1.52 (2785 d. o.f.). 

One of the ionisation phases in lTurner et al.l J2008 ) with log £ 
of 2.2 has a partial covering fraction. Therefore, at this stage of 
the fitting process we tested whether any of our three phases 
is partially covering the source. We freed the covering fraction 
(fcov) of the three phases from the default value of 1 . The cover- 
ing fraction of one of the phases (Phase B with log £ of 2.4) was 
found to change from 1 to 0.55, whist the fcov of the other two 
phases remained practically unchanged. By freeing the covering 
fraction of phase B,xl fell from 1.52 (2785 d.o.f.) to 1.44 (2784 
d.o.f.). In fact there are absorption lines in the 1 1-13 A region of 
the RGS spectrum which are not fitted well if fcov is fixed to 1, 
as well as the UTA and the O vn He-a line. So a partial covering 
of the continuum is necessary for phase B, but this is not the case 
for the other two phases. Therefore, in the following we allowed 
phase B to have a partial covering fraction and fcov of phases A 
and C were fixed to 1 . A plot of the whole RGS spectrum, the 
EPIC-pn data and the best fit model for Obs. 2 is shown in Fig. 




Table 3. Best fit parameters of the continuum, obtained from the three-phase xabs model simultaneous fit to the EPIC-pn, RGS 1 and RGS2 spectra. 



Obs. 


Power-law: 
Photon Index (r) Normalisation^ 


Modified black body: 
Temperature^ Normalisation^ 


1 


1.85 ±0.02 


3.3 ±0.1 


189 ±2 


1.3 ±0.1 


2 


1.82 ±0.01 


2.8 ±0.1 


186 ±3 


1.5 ±0.1 


3 


1.70 ± 0.02 


1.9 ±0.1 


211 ±4 


0.5 ±0.1 


4 


1.84 ±0.02 


2.9 ±0.1 


189 ±2 


1.4 ±0.1 



" 10 51 photons s -1 keV" 1 at 1 keV. 
* eV. 

c Emitting area times square root of the electron density in 10 33 cm 1 



Table 4. Best fit parameters of the three-phase xabs model for all four observations, obtained from the simultaneous fitting of the EPIC-pn, RGS 1 and RGS2 spectra. 







Phase A: 








Phase B: 








Phase C: 






Obs. 




mi 


Flow v c 


RMS v c 


log?" 




Flow v c 


RMS v c 


fcovi 


log?" 




Flow v c 


RMS v c 


^/d.o.f.l 


1 


0.95 ± 0.02 


0.33 ±0.01 


-100 ±40 


60 ± 10 


2.41 ± 0.05 


1.7 ±0.2 


-1500 ± 100 


500 ± 100 


0.64 ± 0.09 


3.00 ± 0.05 


1.2 ±0.2 


-800 ± 300 


300 ± 100 


1.41/2746 


2 


0.96 ± 0.02 


0.34 ± 0.01 


-100 ± 30 


40 ± 10 


2.42 ± 0.04 


2.0 ± 0.3 


-1500 ± 100 


500 ± 100 


0.55 ± 0.04 


2.99 ± 0.03 


1.3 ±0.2 


-1000 ± 100 


400 ± 100 


1.44/2784 


3 


0.87 ± 0.02 


0.43 ± 0.02 


-200 ± 40 


40 ± 10 


2.43 ± 0.03 


3.2 ± 0.4 


-1600 ±200 


400 ± 100 


0.64 ± 0.04 


3.07 ± 0.04 


1.7 ±0.2 


-800 ± 200 


400 ± 100 


1.41/2689 


4 


0.97 ± 0.02 


0.37 ± 0.01 


-200 ± 20 


30 ± 10 


2.39 ± 0.05 


2.0 ±0.1 


-1500 ± 100 


600 ± 100 


0.54 ± 0.04 


2.99 ± 0.02 


1.7 ±0.2 


-1000 ± 100 


500 ± 100 


1.42/2724 



" erg cm s 1 . 

* 10 22 cirT 2 . 

c km s -1 . 

d The covering fraction of Phases A and C is 1.0. 

e Thexl values are further improved by including the modelling of the observed emission features discussed in Sect. 15. 51 



Table 5. The elemental abundances of the three-phase xabs model of Obs. 2 fitted simultaneously to the EPIC-p n, RGS1 and RGS2 spectra. The abundances of the three phases 
were coupled to each other in the fits. All the abundances are relative to the proto-solar model of Lodders (2003). The reference element is H. 



Element Abundance relative to H 
N 1.5 ±0.1 

O 1.4 ±0.1 

Fe 1.8 ±0.1 

Other elements 1 (f) 
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The best fit parameters for Obs. 2 (power-law, mbb, xabs and 
abundances) are listed in Tables[3][4]and[5] In Fig. [6] close-ups of 
the RGS spectrum and the best fit model are presented. Having 
obtained the xabs model best fit for Obs. 2, we used it as the 
starting point for each of the other observations, leaving all the 
parameters quoted so far free, with the exception of the abun- 
dances of the xabs phases which were fixed at the Obs. 2 best fit 
(the abundance values are shown in Table[5]). The best fit param- 
eters of the continuum and ionisation phases of Obs. 1, Obs. 3 
and Obs. 4, modelled with the xabs photoionised absorption, are 
also shown in Tables [3] and [4] Note that the xl values in Table 
[4] are further improved by modelling the observed emission fea- 
tures discussed in Sect. 15.51 In the following subsection (Sect. 
I5.3l l. we describe how the Obs. 2 best fit was applied to the Obs. 
3 spectrum, in order to investigate whether the covering fraction 
of phase B is variable or not. 



5.3. Spectral variability in Obs. 3 

As shown in Table [H phase B in Obs. 2 has a best fit covering 
fraction (fcov) of 0.55. Our phase B is very similar to the par- 
tially covering phase (zone 3) of iTurner et al.l (l2008|) in terms o f 
ionisation parameter. We find log £ ~ 2.4 versus Turner et"alT s 
log £ ~ 2.2. However, the column density of our phase B 
(Nu ~ 2.0 x 10 22 cm" 2 ) i s about a factor of 10 smaller than the 
one in ITurner etal .1 (120081) (N» ~ 20 x 10 22 cm 2 ). In Sect.lEH 
we discuss the values of Nu of warm absorber phases in NGC 
3516 found by different authors. 

ITurner et al.l d20 08) explained the dip in the X-ray lightcurve 
of Obs. 3 by varying only the covering fraction of their zone 3. 
In order to try and reproduce their results, we applied our Obs. 
2 best fit (described in Sect. 15.2b to Obs. 3 by keeping all the 
parameters fixed and freeing only the covering fraction. As for 
Obs. 2, we fitted the EPIC-pn and RGS spectra simultaneously. 
Applying the Obs. 2 best fit to Obs. 3 we obtained an initial^ 2 
of 225 (271 1 d.o.f.). We found that by freeing only the covering 
fraction no good fit was achieved; the covering fraction goes to 
1 and a. xl of 130 (2710 d.o.f.) is obtained. On the other hand, 
when we freed only the continuum parameters, a much better fit 
was obtained, ;^ 2 of 1.76 (2707 d.o.f.). Figures [7] and [8] depict 
the difference between the two approaches to Obs. 3, proving 
that changing the continuum parameters is a much more feasible 
solution than changing only the covering fraction of phase B in 
order to fit the Obs. 3 spectrum. 

We then tried freeing both the continuum and covering frac- 
tion. In this case, the covering fraction becomes 0.70 and^ 2 goes 
to 1.65 (2706 d.o.f.). We then freed the ionisation parameter and 
Nu of the three xabs phases, parameters of the Fe Ka line and 
the relative normalisations of the EPIC-pn and RGS instruments; 
the xl improved to 1.42 (2695 d.o.f.). Finally, we freed the ve- 
locities of the xabs phases, keeping the abundances fixed. In this 
case, the fit does not improve significantly, a xl of 1.41 (2689 
d.o.f.) was obtained. The best fit results for Obs. 3, as for the 
other three observations, are shown in Tables [3] and [4] 

It is worth noting the advantage of modelling both the RGS 
and EPIC-pn together, which leads to simultaneously constrain- 
ing the broadband continuum using the EPIC-pn and modelling 
the absorption features using the RGS. This approach has a 
drawback in that, due to smaller statistical fractional errors as- 
sociated with the EPIC-pn compared to the RGS, the best fit 
favours the EPIC-pn. However, by making sure that the RGS ab- 
sorption features are fitted properly, this weakness is not a prob- 
lem. 




30 31 32 

Observed Wavelength (A) 



Fig. 6. The Obs. 2 best fit xabs model, fitted simultaneously to 
the RGS and EPIC-pn spectra. The RGS data are shown in blue, 
the EPIC-pn data in red and the models in black. 
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C5 ts s [ Wavelength (A) Observed Wavelength (A) 

Fig. 7. The simultaneous fit made to the EPIC-pn and RGS spec- Fig- 8. The simultaneous fit made to the EPIC-pn and RGS spec- 
tra of Obs. 3, by using the best fit model of Obs. 2 and fitting only tra of Obs. 3, using the best fit model of Obs. 2 and fitting only 
the covering fraction of phase B. Thexl is 130 (2710 d.o.f.). The *e continuum parameters (maintaining fcov = 0.55). The xl is 
RGS data are shown in blue, the EPIC-pn data in red and the 1-76 (2707 d.o.f.). The RGS data are shown in blue, the EPIC-pn 
models in black. data in red and the models in black. The significance of the O vn 

He-y and Caxiv lines is discussed in Sect. 15.41 
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Obs. 2 




Observed wavelength (A) 

Fig. 9. Our best fit warm absorber xabs model showing the three 
phases applied separately to the continuum (phase A: black, 
phase B: blue, phase C: red) from simultaneous fitting of the 
EPIC-pn and RGS spectra of Obs. 2 (top) and Obs. 3 (bottom). 
The best fit parameters of the models are given in Tables [3] @] 
and|5] The models for Obs. 1 and Obs. 4 look identical to that of 
Obs. 2 (top). 

5.4. The partially covering phase B 

Figure [9] shows how each one of the three phases in our warm 
absorber model contributes to the best fit for the four observa- 
tions. Phases A, B and C are plotted in black, blue and red, 
respectively. Phase A, which is the lowest ionisation phase, is 
associated with M-shell Fe absorption form ing the Unresolved 
Transition Array (UTA) (iBehar et al .11200 lh between 16 and 17 
A. Phase C, the highest ionisation phase, is associated with 
highly ionised Fe. Phase B, of intermediate ionisation, is the 
phase with the partial covering fraction. 

In the previous subsection (Sect. 15.3b , we showed that by 
changing only the covering fraction of phase B, one cannot 
fit the Obs. 3 spectrum. To confirm that fcov of phase B does 
not change significantly between observations, we examined our 
model closely to look for absorption lines unique to phase B (i.e. 
that do not appear in phases A and C) and to see if they change 
between Obs. 2 and Obs. 3. This is important because by ob- 
serving a change in the depth of these lines, we can investigate 




17 18 19 20 21 22 23 24 25 



Observed wavelength (A) 

Fig. 10. The Obs. 2 best fit xabs model, showing that the O vn 
He-y and Caxiv absorption lines are only modelled in phase B. 
Phase A is shown in black, phase B in blue and phase C in red. 



the effect of a covering fraction change. We found three of such 
absorption lines, which are only modelled by phase B. They are 
O vn He-y (17.8 A), Ca xiv (24. 1 A) and S xn (36.4 A). Figure[TU] 
shows the presence of O vn He-y and Ca xiv lines in our model; 
the S xn line is not displayed separately to save space, but the fol- 
lowing is also true for the S xn line. Comparing the fits in Figs. 
[6] (Obs. 2) and [8] (Obs. 3), we can see that the Ovu He-y and 
Ca xiv are fitted well in both. It is important to note that in both 
figures, the covering fraction value is the same (i.e. 0.55). If the 
covering fraction of phase B had changed significantly in Obs. 3, 
we would have expected the Ovu He-y and Caxiv lines in Fig. 
[8]not to be fully fitted. The fact that they are fitted well by vary- 
ing only the continuum parameters is an extra piece of evidence 
that the variability seen in the lightcurve is mostly intrinsic to 
the source, and not the warm absorber. 

5.5. Narrow and broad emission features 

The only significant narrow emission feature in the RGS spectra 
of the four observations of NGC 3516 is the O vn forbidden line 
(rest frame wavelength of 22.10 A). Figure QT] shows the Obs. 
2 and Obs. 3 spectra and best fit around the O vn triplet region. 
This line was modelled using a Gaussian profile (gaus) and its 
best fit parameters are given in Table [6] The RGS spectrum of 
NGC 35 16 around t he O vn triplet is very sim ilar to that of NGC 
5548 presented bv ISteen brugge et al. (2003): in that case, the 
forbidden line is present in emission and outflowing at around 
1 10 km s _1 , the intercombination line is absent or very weak and 
the resonance line is seen in absorption; there are also absorptio n 
lines by O v and O vi. As suggested bv lSteenbrugge et al] (l2003). 
one possible explanation for not detecting O vn intercombination 
emission is blending of the intercombination line with the O vi 
absorption line. 

From Figs. [6] and [8] an excess of emission is noticeable in the 
RGS data at around 18.0 to 18.5 A. This excess is observed in all 
four observations. Also, in this region of the spectrum, phase A, 
the lowest-ionised phase, produces absorption by O vi. We first 
fitted the excess in Obs. 2 with a simple Gaussian profile; best 
fit line wavelength of 18.16 ± 0.04 A (in the rest frame of the 
AGN) and FWHM of 0.47 ± 0.09 A (cr,, = 3300 + 600 km s" 1 ) 
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Table 6. The best fit parameters of the narrow O vn forbidden emission line shown in Fig. Q~T] obtained by adding a gaus component 
to the best fit xabs model of all four observations. The errors correspond to a Ajf 2 of 2. 



Parameter Obs. 1 Obs. 2 Obs. 3 Obs. 4 

A ^ 22.10 22.10 22.10 22.10 

A^ 22.11 ±0.06 22.09 ±0.04 22.11 ±0.02 22.12 ±0.08 

Flowv^ +100 ± 800 -100 ±500 +100 ±300 +300 ±1100 

FWHM^ <0.1 <0.1 <0.1 <0.1 

cr v ' < 600 < 600 < 600 <600 



Line Normalisation ' 1.0 ±0.6 1.0 ±0.6 1.7 ±0.5 0.9 ± 0.7 



" Theoretical rest frame wavelength in A. 

* Wavelength in the rest frame of NGC 35 16 in A. 

c km s _1 . 

" A. 

e RMS velocity in km s -1 (cr v = FWHM/ Vln 256). 
f 10 49 photons s-'. 




21.5 22 
Observed Wavelength (A) 



Obs. 3 


s 

CJ 
u 

X 

> 
o 






i|tt^. 



21.5 2 
Observed Wavelength (A) 



Fig. 11. The presence of forbidden O vn line emission in Obs. 2 
(top) and Obs. 3 (bottom). The RGS data are shown in blue, the 
EPIC-pn data in red and the models in black. 



were obtained with °f 1-40 (2781 d.o.f.). The wavelength of 
the emission line does not correspond to any prominent transi- 
tion, unless this is emission by O vm Ly-a (rest-frame A of 18.97 
A) with an outflow velocity of about 13000 km s _1 . Therefore, 





18.5 

Observed Wavelength (A) 



Fig. 12. Top panel: the Obs. 2 best fit xabs model fitted simul- 
taneously to the RGS and EPIC-pn spectra with the relativistic 
emission line included (compare with Fig.|6]where the line is not 
included). The RGS data are shown in blue, the EPIC-pn data in 
red and the models in black. Middle panel: residuals of the fit in 
the top panel. Bottom panel: the Obs. 2 best fit model of the rel- 
ativistic O vm Ly-a emission line added to the continuum. The 
model is shown in the observed frame, and for clarity of presen- 
tation no absorption component of our model is included in this 
figure. 



as an alternative we investigated whether the excess could be 
caused by a relativistically-broadenedO vm accretion disk emis- 
sion line. If emission from the inner disk is influenced by general 
relativistic effects, the lines would appear to be broad and asym- 
metric. Such lines were first proposed to explain the spectra of 
the Narrow Line Seyfert 1 galaxie s MCG -6-30-15 and Mrk 766 
dBranduardi-Ravmont et al.ll2.Q0lh . We added a narrow emission 
line (delt component in SPEX) at the rest frame wavelength of 
the O vm Ly-a line (18.97 A) at the redshift of NGC 3516, and 
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convolved it with the relativistic disk line profile of lLaorl (Il99ll) 
(laor). Parameters of the laor component are the disk inner and 
outer radii r; n and r out , its inclination angle i and the emissivity 
index q, which appears in the emissivity law, taken to be pro- 
portional to (r 2 + h 2 y q l 2 , where r is the radius in the accretion 
disk and h a scale height dKaastra et alj|2002l) . By adding the 
relativistic O vm line to the xabs best fit model of Obs. 2, our fit 
improved as^ 2 fell from 1.44 (2784 d.o.f.) to 1.40 (2779 d.o.f.) 
and the excess was well fitted. The line profile and the best fit 
are shown in Fig.[T2l 

The inner disk radius is found to be 14 + 3 GM/c 2 and the 
outer disk radius is poorly constrained. The signal of the red tail 
of the line is not strong enough to be modelled accurately; there- 
fore, as the asymmetry of the line is not very clear in the data, the 
broadened Gaussian and the Laor profile fit the excess equally 
well in our modelling; however, as explained earlier, a Gaussian 
profile does not correspond to any prominent transition. 

The disk inclination angle is constrained very well in our 
model for all four observations; this is because its value very 
much depends on the position of the bulk of the excess emission 
in the spectrum: a higher inclination angle causes the emission 
line profile, shown in Fig.[T2](bottom), to move to lower energies 
and vice versa. As the excess emission is clear with respect to the 
continuum, the disk inclination angle is derived to be 33° + 2. 

6. Discussion 

In this work, we have demonstrated the advantage of combining 
RGS high resolution with EPIC-pn spectroscopy in modelling 
an AGN warm absorber. EPIC-pn constrains the continuum well, 
but has insufficient resolution to resolve absorption and emission 
features seen in the warm absorber (e.g. the Fe UTA and Ovm 
Ly-a absorption as well as the O vn He-a (f) emission line shown 
in Figs. l6l and fTTb. Therefore, to investigate changes in the warm 
absorber between observations, careful examination and fitting 
of the RGS spectrum are needed alongside fitting the EPIC-pn 
spectrum. In the following subsections we discuss our results 
and their physical implications. 

6.1. The FeKaline 

In this work we modelled the broad Fe Ka line at ~ 6.4 keV 
with a simple Gaussian profile. For an in- depth spectra l analy sis 
of the Fe K regime, we refer the reader to lTurner et al.l d2008l) in 
which the Chandra HETG spectrum is used to model the Fe Ka, 
Fe KB and Fe xxvi emission lines in detail. However, from our 
simple modelling of the line, described in Sect. [4] (Table [1), we 
find that the Fe Ka flux remains practically unchanged between 
observations, whereas the power-law flux falls significantly in 
Obs. 3. It is known that the variability of the Fe Ka line in gen- 
eral is no t correlated to that of the observed continuum in a trivial 
manner (Miniutti & Fabi anl2004l) . The iron line does not always 
respond to variations in the continuum: in some cases the Fe Ka 
line can appear to b e constant while the co ntinuum varies by a 
large amplitude (e.g. lMarkowitz et al.l2003l) . Using the EPIC-pn 
spectra, we cannot examine the Fe Ka profile so accurately as 
with the Chandra HETG, so we will not discuss this any further. 

6.2. The warm absorber structure 

We have found that three phases of ionisation are adequate to 
model the warm absorption in NGC 3516. The elemental abun- 
dances are unlikely to be different in the three phases, since it is 



possible to fit the spectrum very well by coupling all the abun- 
dances at values close to solar (Table [5]). There is also no need 
for a large over-abundance of iron to fit the UTA. 

A comparison of A^h and outflow velocity values of warm 
absorber phases from recent observations of NGC 35 16 reported 
by different authors is shown in Table [7] The ionisation parame- 
ter s of our phases B an d C are remarkably similar to those found 
bv lTurner et al.l d2005l) during the 2001 observations. We do not 
find requirement for absorption by a fourth phase ( identified by 
Chan dra HETG to correspond to log £ ~ 4.3 in iTurner et al.l 
2008) from our simultaneous EPIC-pn and RGS fittings, which 
are not very sensitive to such high ionisation gas. Note that log £ 
values are not strictly comparable as they are SED dependent, al- 
though the differences are likely to be only of the order of some 
tenth of dex. The outflow velocity of ou r partially covering phase 
B is practically the same as the one in Turn er et al.l J20 08) (their 
partially covering Zone 3); these phases also have very similar 
ionisation parameters, suggesting t hey represent the same pho- 
toionised gas. The second phase in ITurner et al.l d2005l) . which 
was also partially covering, is again very similar to our phase B 
in terms of ionisation parameter and outflow velocity (albeit this 
is smaller in 2005). The pattern that emerges by comparing the 
outflow velocities of all the phases found by different authors 
is the existence of low-ionised phases with very small outflow 
velocity and higher ionisation phases with outflow velocities be- 
tween 1000 km s _1 and 1500 km s _1 . 

The partially covering phase is the phase with the highest 
outflow velocity and most likely closest to the central engine. 
Since we do not find changes in the occultation of the nuclear 
source, there is no evidence of a transverse component of veloc- 
ity for the phase. The minimum distance of this phase from the 
central engine can be estimated assuming the outflow velocity 
we measure to be greater than or equal to the escape velocity 
v eS c = y/2GM/r, where G is the gravitational constant, M the 
black hole mass and r the distance of the absorber phase from 
the black hole. Since the partially covering phase B is outflowing 
at 1500 km s _1 , u sing the NGC 3516 b lack hole mass estimate 
of 2.95 x 10 7 M (iNikotaiuk et alj|2006l) . we obtain r > 0.1 pc. 
Using the expression for the ionisation parameter definition (Eq. 
[TJ and taking L; on 3.85 x 10 43 erg s (derived from the fit of 
Obs. 2) we find the hydrogen number density n < 1.3x 10 8 cm -3 . 
Combining this with the relation A^h ~ nhr, where Ar is the 
thickness of the absorber phase, we find Ar > 5 x 10~ 5 pc. This 
suggests the presence of a thin spherical shell of gas, which is a 
feasible scenario for partial covering. 

The column densities A^h of all three phases of our warm ab- 
sorber cover a relatively small range: ~ 0.4-2 xlO 22 crrT 2 , and 
w e do not find "heavy" (> 2 x 10 22 crrT 2 ) absorber phases as 
in ITurner et al.l (120051 120 08), especially those corresponding to 
our phase B. However, our range of A^h is consistent with the 
values obtained from o ther X-ray observations of NGC 3516: 
Costantini et al. | d2000T) found a warm absorber column of ~ 
2-3 XlO 22 cm' 2 from BeppoSAX spectra; iNetzer etal] (120021) 
found a line-of-sight absorber with A^h of ~ 1 x 10 22 crrT 2 from 
both a 1994 ASCA obse rvation and a 2000 Chandra observation; 
Mar kowitz et al.l d20 08) found their two warm absorber phases to 
have A^h of ~ 4 and 5.5 XlO 22 cmT 2 . 

ITurner et al.l (2005) concluded that their three absorber 
phases cannot be in thermal equilibrium with each other in 
NGC 3516. We can con sider the photoion isation equilibrium 
curve technique (S-curve, Krolik et al. 1981) for th e SED used in 
our analysis (NGC 5548, lsteenbrugge et alj|2005b . The S-curve 
marks the points of temperature T versus pressure form of the 
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Table 7. Comparison of the ionisation parameters, hydrogen column densities and outflow velocities of the NGC 3516 warm 
absorber phases found in recent observations, shown in increasing order of 
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-1100 '• 
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26 


-1000 


4 


3.0 
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-900 


C 



" From Apr. and Nov. 2001 XMM-Newton EPIC-pn spectra. The values of the ionisation parameter in £ form are given in lTurner et al.l ( 2008). 

b From Oct. 2005 Suzaku XIS and HXD spectra. 

c From Oct. 2006 XMM-Newton EPIC-pn and Chandra HETG spectra. 

d From Oct. 2006 XMM-Newton EPIC-pn and RGS spectra. 

e erg cm s" 1 . 

/ 10 22 cm" 2 . 

s Flow velocity in km s -1 . 

* Based on the velocity obtained in lTurner et al.l d2005l) . 



ionisation parameter S at which the pressure P is constant. The 
parameter H is related to £, by 



= 0.961 x lO 4 ^ 
4ncr 2 P T 



(2) 



where L; on is the luminosity of the ionising source, c is the 
speed of light, r is the distance of the warm absorber phase 
from the ionising source, P i s the pressure and T the temper- 
ature. [SteenbmggEetaf] d2005h found that the S-curve for NGC 
5548 does not have negative slope anywhere; applying this to 
NGC 35 16, we can conclude that the three absorption phases we 
identify in NGC 3516 outflowing with different velocities are 
unlikely to be in pressure equilibrium. However, it is worth not- 
ing that different photoionisation codes, variations in the shape 
of the SED and selec ting differen t en ergy ranges for the con- 
struction of the SED dSmifhll2009l and ICnakravortv et alj r2009) 
can produce some differences in the S-curve. Because the ion- 
isation parameters, column densities and outflow velocities of 
our phases B and C are similar (although not in pressure equi- 
librium), and more importantly the covering fraction of phase B 
is not variable, phase B could have the same origin as the fully 
covering phase C. This disfavours a clumpy disk wind scenario 
as the most likely explanation for the origin of phase B. 

6.3. Narrow and broad emission features 

The RGS spectrum shows evidence for both narrow and broad 
emission features. The Ovn (f) emission line is found to have a 
lower flow velocity compared to the observed absorption lines 
of the high ionisation phases B and C; the flow velocity is in 
fact consi stent with zero, an d with that of the lowest ionisation 
phase A. Smith et al] (120081) have observed the O vn (f) line to 
have a very low outflow velocity compared to other soft X-ray 
emission lines in Ark 564, although in that case all absorption 
phases also have flow velocities consistent with zero. This lack 
of a velocity shift from the rest wavelength is in agreement with 
a scenario in which we are observing emission from a fully vis- 
ible spherically symmetric outflowing shell around the nuclear 
source, thus no net velocity of O vn will be detected. Another 
possible explanation is that the O vn line originates in the Narrow 
Line Region (NLR), where the emission lines are produced in 
a gas of compa ratively small velocity dispersion and low den- 
sity (see for e.g. Kaastra et al.l d2002l) . where the forbidden O vn 



line in NGC 5548 is suggested to come from the NLR). As the 
forbidden line is the most intense line in the O vn triplet, a pho- 
toionised gas with low density is required. 

We have found possible evidence for a relativistically broad- 
ened emission line of O vra Ly-a. However, the line is weaker 
tha n the strong Ovm, Nvn and C vi Ly-a lines repor ted 
bv iBranduardi-Ravmont et al l (1200 ll) and ISako et ail d2003l) in 
MCG -6-30-15 and Mrk 766. Evidence for weak broadened 
emission l ines of O vi ii and N vn Ly-a is also found in NGC 
5548 by Kaast ra et al.l d2002l) . The disk inclination angle ob- 
tained f rom our fit to the O v m Ly-a profile is 33°±2. IWu & Hani 
(l200ll) : IZhang & Wul d2002l) have calculated the inclination an- 
gles of the Broad Line Region (BLR) of several Seyfert 1 galax- 
ies using already published bulge stellar velocity dispersions and 
black hole masses estimated by reverberation mapping. They 
calculated an inclination angle of 38.3° ± 7.6 for NGC 3516, 
which would imply co-planarity of the BLR with the AGN ac- 
cretion disk. 



6.4. Intrinsic continuum versus warm absorber variability 

NGC 35 16 has a histor y of large amplitude continuum variability 
between observations. iNetzer et al.1 d2002l) reported a large drop 
in flux (factor of ~ 50 at 1 keV) between an ASCA observation 
in 1994 and a Chandra observation in 2000. They found that 
the observed flux and spectral variability at these epochs were 
consistent with a constant column density of line-of-sight mate- 
rial reacting to changes in the ionising continuum. In the 2006 
XMM-Newton data we have found the observed spectral and 
flux variability to be unrelated to changes i n the covering frac - 
tion of phase B, unlike what is reported bv lTurner et al.l (12008). 
Furthermore, the X-ray absorption line depths are sensitive to 
changes in the covering fraction; from close examination of RGS 
high resolution spectra, there is no evidence to suggest that the 
covering fraction of phase B changed between observations. 

The only parameters which indicate a change during the 
low flux XMM-Newton Obs. 3 are the ionisation parameter £ 
of phase A and the column density Nn of phase B. We find 
log £ ~ 0.9 and Nu ~ 3 x 10 22 cirT 2 in Obs. 3, whereas in 
the other three observations they are ~ 1.0 and ~ 2 xlO 22 cm -2 , 
respectively. However, these changes in the warm absorber pa- 
rameters are very small compared to those in the continuum pa- 
rameters: power-law slope T from ~ 1.8 to 1.7 in Obs. 3, and 
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normalisation from ~ 3 to 2 x 10 51 photons s _1 keV at 1 keV; 
modified black body temperature from ~ 190 to 210 eV in Obs. 
3, and normalisation (emitting area times square root of electron 
density) from ~ 1.5 to 0.5 x 10 33 cm - 5 . 

The power-law contribution to the continuum heavily out- 
weighs the modified black body in all our four observations and 
is responsible for nearly all of the continuum variability. From 
the X-ray lightcurve (Fig. [1} and spectrum (Fig. the varia- 
tion in Obs. 3 is larger in the soft X-ray (0.2-2.0 keV) energy 
band than in the hard X-ray (2.0-10.0 keV) band. This type of 
variability has been seen in ot her X-ray observat ions of Seyfert 
1 AGN (such as NGC 7469, iBlustin et alj|2003b in which the 
source is softer when brighter. 

7. Conclusions 

1. We have studied the warm absorber in NGC 3516 by 
analysing in detail EPIC-pn and RGS high resolution spec- 
tra from four observations made in October 2006 by XMM- 
Newton. The warm absorber consists of three phases of ion- 
isation: phase A (log £ ~ 0.9), phase B (log £ ~ 2.4) and 
phase C (log £ ~ 3 .0) in increasing order of ionisation. Phase 
A has a hydrogen column density of ~ 0.4 x 10 22 cm -2 , 
which is smaller than those of the other two phases (~ 
1 -2x 10 22 cm -2 ). There is evidence that the lower-ionisation 
phase A is outflowing at ~ 100 kms -1 , whereas the two 
higher-ionisation phases are outflowing faster at around 1000 
to 1500 km s" 1 . 

2. Phase B covers about 60% of the source continuum. We in- 
vestigated whether variation in the covering fraction of phase 
B could account for the o bserved flux and sp ectral variability 
in Obs. 3 (as claimed bv lTurner et al.ll2008l) . We found that: 
(1) the covering fraction does not show significant variation 
between observations; (2) even if the covering fraction is al- 
tered significantly, this does not properly account for the ob- 
served variability in the EPIC-pn and especially RGS spec- 
tra. This makes a clumpy disk wind scenario a rather unfea- 
sible expla nation. 

3. Similar to Netzeretal ] (12002b . we conclude that the vari- 
ably in the 2006 observations presented here (albeit much 
smaller than in previous observations) is better understood as 
the consequence of changes in the source continuum emis- 
sion than in the warm absorber. Our results suggest that the 
X-ray variability of NGC 3516 (and by inference, possi- 
bly that of other AGN where a similar behaviour has been 
observed) cannot be reduced to occultation and absorption 
effects, su ch as proposed bv ITurner et al.l (120081) (see also 
review by ITurner & Milled [2009, and references therein); 
rather, the variability is likely to arise in a scenario where 
intrinsic changes of the continuum and ionisation state are 
also important. In this context, a careful analysis of the soft 
X-ray spectrum at high resolution, in combination with data 
in the Fe Ka regime, can provide essential constraints and 
clues to the source physical behaviour. 
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